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Abstract. The Riccati equations reducible to first-order linear equations 
by an appropriate change the dependent variable are singled out. All these 
equations are integrable by quadrature. 
^ ', A wide class of linear ordinary differential equations reducible to al- 

gebraic equations is found. It depends on two arbitrary functions. The 
method for solving all these equations is given. The new class contains the 
constant coefficient equations and Euler's equations as particular cases. 
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On 
OO 

The old-fashioned title of the present paper indicates that it is dedi- 
cated to quite old topics. Namely, it deals with a problem on integration by 
quadrature of Riccati equations investigated, in terms of elementary func- 
tions, by Francesco Riccati and Daniel Bernoulli some 280 years ago for the 
special Riccati equations (see, e.g. [1]) 
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y' = ay 2 + bx a , a, b, a = const., 

H [ 

and with an integration of higher-order linear equations by reducing them 
to algebraic equations. The later property was discovered by Leonard Euler 
in the 1740s for the constant coefficient equations 

y (n) + A iy {n - l) + ■■■ + A n - iy ' + A n y = 0, A 1} ...,An = const., 

as well as for the equations of the form 



n (n) 



x y 



+ A^-y™" 1 ) + • • • + A n ^xy' + A n y = 0, A 1 ,...,A n = const., 



known as Euler's equations. 

It will be shown in what follows that these classical results can be ex- 
tended to wide classes of equations. 



^irst published in Archives of ALGA, vol. 5, 2008, pp. 27-62. 
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Chapter 1 
Riccati equation 



1 Introduction 

Consider the special Riccati equation 

y' = ay 2 + bx a , a, b, a = const. (1.1) 
If a = 0, Eq. (1.1) is integrable by the separation of variables: 

dy 



ay 2 + b 



= dx. 



Another easily integrable case is a = — 2. Then the change of the dependent 
variable 

1 

z — — 

y 

maps Eq. (1.1) to the homogeneous equation 

dx . (x) . 

which is integrable by quadrature. 

F. Riccati and D. Bernoulli noted independently that Eq. (1.1) can be 
transformed to the case a — 0, and hence integrable by quadrature in terms 
of elementary functions if a takes the values from the following two series: 
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The series (1.2) are given by the formula 

4k 

a = with fc = ±l, ±2, .... (1.3) 

2k -1 

It is manifest from (1.3) that both series in (1.2) have the limit a = —2. 
For a derivation of the transformations mapping Eq. (1.1) with a having 
the form (1.3) to an integrable form, see [1], Chapter 1, §6. 

J. Liouville showed in 1841 that the solution to the special Riccati equa- 
tion (1.1) is integrable by quadrature in terms of elementary functions only 
if a has the form (1.3). 

It is well known that the general Riccati equation 

y' = P(x) + Q(x)y + R(x)y 2 

can be rewritten as a linear second-order equation. But this kind of lin- 
earization by raising the order does not solve the integration problem. 
Therefore, I will investigate a possibility of linearization of the Riccati equa- 
tions without raising the order and will show that all Riccati equations of 
this type can be integrated by quadrature. 

2 The linearizable Riccati equations 

The following theorem is closely related to the theory of nonlinear superpo- 
sitions discussed in [2], Section 6.7. 

Theorem 2.1. The first-order ordinary differential equation 

y' = f(x,y) (2.1) 
can be reduced to a linear first-order equation 

dz 

— =p{x) + q{x)z (2.2) 
ax 

by a change of the dependent variable y, 

z = z(y), (2.3) 
if and only if Eq. (2.1) can be written in the form 

y 1 = T l {x)Uy)+T 2 {x)Uy) (2.4) 

such that the operators 

* = m^ x * = Uv)^ (2 - 5) 

span a two-dimensional (or a one-dimensional if X\ and X 2 are linearly de- 
pendent) Lie algebra called in [3] the VGL (Vessiot-Guldberg-Lie) algebra. 
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Proof. Let Eq. (2.1) be linearizable. Then we can assume that it is already 
reduced by a certain change of the dependent variable 



z = C(y) (2.6) 

to a linear equation (2.2), 

^ = p(x) + q(x)z. 
ax 

The VGL algebra associated with Eq. (2.2) is two-dimensional and is 
spanned by the operators 

x °=4- (2 - 7) 

The form of Eq. (2.4) and the algebra property [Xl, X2] = aX\ + PX 2 
remain unaltered under any change (2.3) of the dependent variable. There- 
fore, rewriting Eq. (2.2) in the original variable y = C _1 (z) obtained by the 
inverse transformation to (2.6), we arrive at an equation of the form (2.4) 
for which the operators (2.5) span a two-dimensional Lie algebra. Since the 
equation obtained from Eq. (2.2) by the inverse transformation to (2.6) is 
the original equation (2.1) we have proved the "only if part of the theorem. 

Let us prove now the "if part of the theorem. Namely, we have to 
demonstrate that any equation of the form (2.4) such that the operators 
(2.5) span a two-dimensional Lie algebra, is linearizable. If the operators 
(2.5) are linearly dependent, then C, 2 { x ) — 7 — const., and hence 

Eq. (2.4) has the form 

y' =[T l {x)+ 1 T 2 {x)]i l {y)- 
It can be reduced to the linear equation 

z' = T 1 (x)+ 1 T 2 (x). 

upon introducing a canonical variable z for 

by solving the equation X\(z) = 1. 

Suppose now that the operators (2.5) are linearly independent. It is 
clear that Eq. (2.9) will be linearized if one transforms the operators (2.5) 
to the form (2.7). One can assume that the first operator (2.5) has been 
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already written, in a proper variable z, in the form of the first operator Xi 
given in (2.7): 

az 

Let the second operator (2.5) be written in the variable z as follows: 

*.=/«!■ 

We have 

[X l5 X 2 ] = /'(^^ 

and the requirement [Xi,X 2 ] = aX\ + /3X 2 that X±,X 2 span a Lie algebra 
L 2 yields the differential equation 

f' = a + l3f, 

where not both a and j3 vanish because otherwise the operators X\ and X 2 
will be linearly dependent, hence f'(x) 7^ 0. Solving the above differential 
equation, we obtain 

/ = az + C if = 0, 
f = C^ z -% if (3^0. 

f = ax + C X 2 = oa:-^- + CXi, if 6 = 0, 

dx 

/ = Ce ^_a ^ x 2 = Ce hx ^--yX 1 , if 6^0. 
o dx b 

In the first case we have 

d 

X 2 = az— + CX U 

az 

and hence a basis of L 2 is provided by (2.7). In the second case we have 

and one can take basis operators, by assigning /3z as new z, in the form 

X - 9 X - e z 9 

Finally, substituting z = e~ z we arrive at the basis (2.7), thus completing 
the proof of the theorem. 
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The following theorem characterizes all Riccati equations that can be 
reduced to first-order linear equations by changing the dependent variable 
(see [4], Russian ed., Theorem 4.3; [3], Section 11.2.5 and Note [11.4]; see 
also Theorem 3.2.2 in [2]). 

Theorem 2.2. The following two conditions 1° and 2° are equivalent and 
provide the necessary and sufficient conditions for the Riccati equation 

y' = P(x) + Q{x)y + R(x)y 2 (2.8) 

to be linearizable by a change of the dependent variable (2.3), z = z(y). 
1°. Eq. (2.8) has a constant solution y = c including c = oo. 
2°. Eq. (2.8) has either the form 

y' = Q(x)y + R(x)y 2 (2.9) 

with any functions Q(x) and R(x), or the form 

y' = P{x) + Q{x)y + k[Q(x) - kP{x)]y 2 (2.10) 

with any functions P(x), Q(x) and any constant k. 

Proof. The conditions 1° and 2° are equivalent. Indeed, let Eq. (2.8) have 
a constant solution y = c. Then 

= P( x ) + Q(x)c + R(x)c 2 . (2.11) 

If c = 0, Eq. (2.11) yields P(x) = 0, and hence Eq. (2.8) has the form (2.9). 
If c^O, Eq. (2.11) yields 

R(x) = -\ \P(x) + cQ(x)] = -- \Q(x) + - P(x)] . 
c c c 

Denoting k — — 1/c we get 

R(x) = k[Q(x) - kP(x)]. 

Hence Eq. (2.8) has the form (2.10). Thus, we have proved that 1° => 2°. 
Conversely, let Eq. (2.8) satisfy the condition 2°. It is manifest that Eq. 

(2.9) has the constant solution y = 0. Furthermore, one can verify that Eq. 

(2.10) has the constant solution y = —1/k. This proves that 2° ^> 1°. 

Let us turn to the necessary and sufficient conditions for linearization. 
The operators 

d „ , d 



associated with Eq. (2.9) have the commutator [Xi,X 2 ] = X 2 , and hence 
span a two-dimensional Lie algebra. Furthermore, the operators 

X 1 = (1-A;V)|, X 2 = (y + ky^ 

associated with Eq. (2.10) have the commutator [Xi,X 2 ] — Xi + 2kX 2 . 
Hence, they also span a two-dimensional Lie algebra. Therefore, according 
to Theorem 2.1, the condition 2° is sufficient for linearization. 

Note that the equation y' = P{x) + Q(x)y can be regarded as a particu- 
lar case of Eq. (2.10) with k — 0. Since Eq. (2.10) has the constant solution 
y — — 1/k, we conclude that the linear equation y 1 = P(x) + Q(x)y has the 
constant solution y = oo. We conclude that any linearizable equation has 
a constant solution because the change of the dependent variable z — z(y) 
maps a constant solution into a constant solution of the transformed equa- 
tion. Hence, the condition 1° is necessary for linearization. This completes 
the proof of the theorem due to the equivalence of the conditions 1° and 2°. 

Now we will use Theorem 2.2 for integrating the linearizable Riccati 
equations. We will find linearizing transformations for the equations (2.9) 
and (2.10). We will assume that k in Eq. (2.10) is a real number. 



3 Linearization and integration of Eq. (2.9) 

Invoking Eqs. (6.7.5), (6.7.6) from [2], replacing there t and x % by x and 
y, respectively, and identifying Ti(t) and T 2 (t) with Q(x) and R(x), respec- 
tively, we see that the VGL (Vessiot-Guldberg-Lie) algebra associated with 
Eq. (2.9) is a two-dimensional algebra L 2 spanned by the operators 

d 9 d 

Their commutator is [Xi,X 2 ] = X 2 . Introducing the new basis 

X[ = X 2 , X 2 = —X\ 

i.e. taking 

we have a basis in L 2 satisfying the commutator relation 

[X[,X' 2 ]=X' 1 . (3.2) 
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Let us find a change of the dependent variable, z = z(y), such that Eq. 
(2.9) becomes a linear equation (2.2), 

^ = p(x) + q(x)z. 
ax 

The VGL algebra of this equation is spanned by the operators (2.7), 

— d — d 

X\ — — , X 2 — z— 

oz oz 

whose commutator has the same form as Eq. (3.2), i.e. [Xi,X 2 ] = X 1 . 
Consequently, the linearizing transformation z = z(y) is determined by the 
equations X[(z) = 1, X' 2 (z) = z, or 



,dz 



dz 



ay ay 
Integrating the first equation (3.3) we obtain 



(3.3) 



z = h A. 

y 

Substituting in the second equation (3.3) we get A = 0. Thus, the linearizing 
transformation is 

, = -i. (3.4) 
In this variable, Eq. (2.9) becomes the linear equation 

z ' = R(x) - Q(x)z, (3.5) 

whence 

Z =[C + J R{x)e$ Q{x)dx dx\ e --l' Q{x)dx , C = const. 
Substituting in Eq. (3.4) we finally arrive at the solution to Eq. (2.9): 



y = - 



C+ R(x)ef Q(x)dx dx 



^fQ(x)dx 



(3.6) 



Remark 3.1. In Section 12.1 the solution (3.6) is obtained by an alterna- 
tive method. 
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Example 3.1. Consider the equation 



y' = - + R(x)y 2 

x 



By Eq. (3.6) yields its solution 



x 

y = 



C — f xR(x)dx 
Example 3.2. The equation 

y' = - + 3xy 2 
x 

is a particular case of the previous equation. Its solution is 

x 



y 



C-x 2 

Example 3.3. The equation 

1 2 

y =y+ -y 

X 



has the solution 



X 

e*^ f e^ 

y = — — — , , , where Ei(x) = — / — dt. 
y C + Ei(x) ' J t 



4 Linearization and integration of Eq. (2.10) 

Now we identify the coefficients P(x) and Q(x) of Eq. (2.10), 

y' = P{x) + Q{x)y + k[Q(x) - kP(x)]y\ (2.10) 

with the coefficients Ti(t) and ^(t) of Eqs. (6.7.5), (6.7.6) from [2] and 
associated with Eq. (2.10) the VGL algebra spanned by the operators 

X 1 = (1-A;V)|, X 2 = (y + ky^. 

Their commutator is [Xl, X 2 ] — Xi + 2kX 2 . We assume that k ^ 0, since 
otherwise Eq. (2.10) is already linear. Therefore, setting 

X[ = X l + 2kX 2 , X' 2 = -X l /(2k) 
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we obtain the new basis 



(4.1) 



satisfying the commutator relation (3.2), [JT(, JT^] = X[. 

The transformation z = z(y) of the operators (4.1) to the form (2.7) is 
determined by the equations 

X[(z) = l, X' 2 (z)=z, 



or 

Integrating the first equation (4.2) we obtain 

, = A 1 



z. 



k(l + ky) 



(4.2) 



Substituting in the second equation (4.2) we get A = l/(2k). Thus, the 
linearizing transformation is 



Eq. (4.3) yields: 



y 



ky - 1 
2k(ky + l) 

l + 2kz 



k(l-2kz) 

Substituting (4.4) in Eq. (2.10) we arrive at the linear equation 



(4.3) 



(4.4) 



z' = —Q(x)+[Q(x)-2kP(x)]z, 



whence 



2k 



C+ J Q ix)e fl2kP(*)-Q(x)]d* dx 



J[Q{x)-2kP(x)]dx 



(4.5) 



(4.6) 



Substituting (4.6) in (4.4) we will obtain the solution to Eq. (2.10). 

Remark 4.1. We have assumed that Q(x) ^ because otherwise Eq. 
(2.10) is separable, 

y' = P(x)(l-k 2 y 2 ). 
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Remark 4.2. Eq. (2.10) has the constant solution 

y = -\- (4.7) 

We excluded the singular solution (4.7) in the above calculations and as- 
sumed that 1 + ky ^ 0. 

Example 4.1. Let us integrate the equation 

y' = x + 2xy + xy 2 . (4.8) 

Here 

P{x) = x, Q(x) =2x, k = 1. 

Therefore the linearized equation (4.5) is written z' = x. Integrating it and 
denoting the constant of integration by C/2, we obtain 

z = \{C + x 2 ) 

in accordance with Eq. (4.6). Substituting the above z in Eq. (4.4) we 
obtain the following solution to Eq. (4.8): 



^ + C + x 2 , An , 



Remark 4.2 gives also the singular solution y — — 1. 

Remark 4.3. Eq. (4.8) can also be integrated by separating the variables. 
Example 4.2. The equation 

y ' = x 2 + {x + x 2 )y + ^(2x + x 2 )y 2 (4.10) 

has the form (2.10) with 

P(x) = x 2 , Q{x) = x + x 2 , k — - ■ 

The linearized equation (4.5) is written 

z' = x + x 2 + xz 

and has the solution 



(c+ J{x + x 2 )e' x2 ' 2 dx^ e x2 / 2 . 



11 



Substitution in Eq. (4.4) yields 

y = 2- 

1 — z 

Hence, the solution to Eq. (4.10) is given by 



1+(C + J(x + x 2 )e~ x2 / 2 dx^j e* 2 /2 
V = 2 1 - (C + J(x + a; 2 )e- 2 /2 dx ) e- 2 / 2 ' (4U) 
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Chapter 2 



Higher-order linear equations 
reducible to algebraic equation 

5 Introduction 

The linear ordinary differential equations with constant coefficients 

j/W + V n - 1 ' + - + A,.i3/' + A n y = 1 Ai,..., A n = const., (5.1) 

and Euler's equations 

+ A lX n - l y^ + ■■■ + A n _ lX y' + A n y = 0, A t , . . . , A n = const., 

(5.2) 

are discussed practically in all textbooks on differential equations. They are 
useful in applications. The most remarkable property is that both equations 
(5.1) and (5.2) are reducible to algebraic equations. Namely, their funda- 
mental systems of solutions, and hence the general solutions can be obtained 
by solving algebraic equations. Then one can integrate the corresponding 
non-homogeneous linear equations by using the method of variation of pa- 
rameters. It is significant to understand the nature of the reducibility and 
to extend the class of linear equations reducible to algebraic equations. 

In the present paper we will investigate this problem and find a wide 
class of linear ordinary differential equations that are reducible to algebraic 
equations. The new class depends on two arbitrary functions of x and 
contains the equations (5.1) and (5.2) as particular cases. The method for 
solving all these equations is given in Sections 8 and 10, and illustrated by 
examples in Section 9. 

The main statement for the second-order equations is as follows (Section 

8). 
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Theorem 5.1. The linear second-order equations 

P{x)y" + Q(x)y' + R(x)y = F(x) 

whose general solution can be obtained by solving algebraic equations and 
by quadratures, have the form 

0V+ {A + cj)' -2or)(f)y' + (B-Aa + a 2 -0a')y = F(x), (5.3) 

where = (p(x), a = a(x) and F{x) are arbitrary (smooth) functions, and 
A,B = const. The homogeneous equation (5.3), 

<p 2 y" + (A + 0' - 2a)<py' + (B - Aa + a 2 - <pa')y = 0, (5.4) 

has the solutions of the form 

y = e J ^) , (5.5) 

where A satisfies the characteristic equation 

A 2 + AA + 5 = 0. (5.6) 

If the characteristic equation (5.6) has distinct real roots Ai ^ A 2 , the 
general solution to Eq. (5.4) is given by 

y{x) = K l e J *(») +K 2 e J *(«o , K 1 ,K 2 = const. (5.7) 

In the case of complex roots, Ai = 7 + 16, \ 2 = 7 — i0, the general solution 
to Eq. (5.4) is given by 

e J ^T dz . (5.8) 

If the characteristic equation (5.6) has equal roots Ai = A2, the general 
solution to Eq. (5.4) is given by 

e J ^(^ dx , K 1 ,K 2 = const. (5.9) 

The general solution of the non-homogeneous equation (5.3) can be obtained 
by the method of variation of parameters. 

Remark 5.1. The equations with constant coefficients and Euler's equa- 
tion are the simplest representatives of Eqs. (5.4). Namely, setting <p(x) = 
1, a(x) = we obtain the second-order equation with constant coefficients 

y" + Ay' + By = 0, 



y(x) 



K x cos [9 




+ K 2 sin ( 9 




K 1 + K 2 



dx 

W). 
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and Eq. (5.5) yields the well-known formula 

V = e Xx , 

where A is determined by the characteristic equation (5.5). 

If we set (f)(x) = x, <j(x) = 0, we obtain the second-order Euler's equa- 
tion (5.2) written in the form 

x 2 y" + (A + l)xy' + By = 0. 

Then Eq. (5.5) yields the particular solutions for Euler's equations: 

y = x\ 

where A is determined again by the characteristic equation (5.5). For details, 
see Section 9. For other functions 4>(x) and a(x) = 0, Eqs. (5.5) are new. 

6 Constant coefficient and Euler's equations 
from the group standpoint 

For the sake of simplicity, we will consider in this section second-order equa- 
tions 

y" + f(x)y' + g(x)y = 0. (6.1) 

6.1 Equations with constant coefficients 

Let us begin with the equations with constant coefficients 

y" + Ay' + By = 0, A,B = const. (6.2) 

Eq. (6.2) is invariant under the one-parameter groups of translations in 
x and dilations in y, since it does not involve the independent variable x 
explicitly (the coefficients A and B are constant) and is homogeneous in the 
dependent variable y. In other words, Eq. (6.2) admits the generators 

Xl = l' x * = % (6 - 3) 

of the translations in x and dilations in y. We use them as follows [3]. 
Let us find the invariant solution for X = X\ + AX 2 , i.e. 

X = — + Xy— A = const. (6.4) 
ox ay 
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The characteristic equation 

dy dx 

of the equation 

for the invariants J(x, y) yields one functionally independent invariant 



J = ye 



-\x 



According to the general theory, the invariant solution is given by J = 
C with an arbitrary constant C. Thus, the general form of the invariant 
solutions for the operator (6.4) is 

y = Ce Xx , C = const. 

Since Eq. (6.2) is homogeneous one can set C = 1 and obtain Euler's 
substitution: 

y = e Xx . (6.5) 

As well known, the substitution (6.5) reduces Eq. (6.2) to the quadratic 
equation (characteristic equation) 

A 2 + AX + B = 0. (6.6) 

If Eq. (6.6) has two distinct roots, Ai ^ A2, then Eq. (6.5) provides two 
linearly independent solutions 

Vi = e XlX , y 2 = e X2X , 

and hence, a fundamental set of solutions. If the roots are real, the general 
solution to Eq. (6.2) is 

y(x) = K ie XlX + K 2 e X2X , K U K 2 = const. (6.7) 

If the roots are complex, Ai = 7 + 16, A 2 = 7 — i6, the general solution 
to Eq. (6.2) is given by 

y(x) = [Ki cos(ex) + K 2 sm(6x)} e 7X , K U K 2 = const. (6.8) 

In the case of equal roots Ai = A 2 , standard texts in differential equations 
make a guess, without motivation, that the general solution has the form 

y(x) = (Ki + K 2 x) e Xlx , K 1 ,K 2 = const. (6.9) 

The motivation is given in [3], Section 13.2.2, and states the following. 
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Lemma 6.1. Eq. (6.2) can be mapped to the equation z" = by a linear 
change of the dependent variable 

y = a{x)z, a(x)^0, (6.10) 

if and only if the characteristic equation (6.6) has equal roots. Specifically, 
if Ai = A2, the substitution 

y = ze XlX (6.11) 

carries Eq. (6.2) to the equation z" = 0. Substituting in (6.11) the solution 
z = K\ + K 2 x of the equation z" = 0, we obtain the general solution (6.9) 
to Eq. (6.2) whose coefficients satisfy the condition of equal roots for Eq. 
(6.6): 

A 2 -45 = 0. (6.12) 

Proof. The reckoning shows that after the substitution (6.10) Eq. (6.1) 
becomes (see, e.g. [2], Section 3.3.2) 

+ I(x)z = 0, 

where 

I{x)=g{x)-\f{x)- l -f\x) 
and the function cr{x) in the transformation (6.10) has the form 

a(x) = e -y fix)dx . (6.13) 

Hence, Eq. (6.1) is carried into equation z" = if and only if I(x) = 0, i.e. 

f\x)+2f'(x)-4g(x)=0. (6.14) 

In the case of Eq. (6.2), the condition (6.14) is identical with Eq. (6.12), 
and the function a(x) given by Eq. (6.13) becomes 

a(x)=e-^ x . (6.15) 

Furthermore, under the condition Eq. (6.12), the repeated root of Eq. (6.6) 
is Ai = —A/2. Therefore Eq. (6.15) can be written 

a(x) =e Xlx 

and we arrive at the substitution (6.11), and hence at the solution (6.9), 
thus proving the lemma. 
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6.2 Euler's equation 

Consider Euler's equation 

x 2 y" + Axy' + By = 0, A,B = const. (6.16) 

It is double homogeneous (see [2], Section 6.6.1), i.e. admits the dilation 
groups in x and in y with the generators 

x ^ x h A ' 2 = 4' (6 - 17) 

We proceed as in Section 6.1 and find the invariant solutions for the linear 
combination X = Xi + XX 2 : 

d d 

X = x— — h Ay— , A = const. (6.18) 
ox ay 



The characteristic equation 



dy dx 



y x 

of the equation X{J) = for the invariants J(x,y) yields the invariant 

J = y x~ x 

for the operator (6.18). The invariant solutions are given by J = C, whence 

y = Cx x , C = const. 

Due to the homogeneity of Eq. (6.2) we can set C = 1 and obtain 

y = x x . (6.19) 

Differentiating and multiplying by x, we have: 

xy' = Xx x , x 2 y" = A(A — l)x x . 

Substituting in Eq. (6.18) and dividing by the common factor Cx x we 
obtain the following characteristic equation for Euler's equation (6.16): 

A 2 + (A- 1) X + B = 0. (6.20) 

Remark 6.1. According to Eqs. (6.6), (6.20), the characteristic equation 
for Euler's equation written in the form 

x 2 y"+ (A + l)xy' + By = (6.21) 

is identical with the characteristic equation (6.6) for Eq. (6.2) with constant 
coefficients. 
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= 0, 

(6.1) 



7 New examples of reducible equations 
7.1 First example 

Let us find the linear second-order equations (6.1), 

y" + f(x)y' + g(x)y = 0, (6.1) 

admitting the operator 

(7.!) 

where a is any real-valued parameter. Taking the second prolongation of 

Xl = x^ - a^-V|j - [«(« " ^"V + 2*°" Y] A , 
we write the invariance condition of Eq. (6.1), 

X 1 (y" + f(x)y' + g(x)y) 
and obtain: 

x a ~ 2 [x 2 f + axf - a(a - l)]y' + x a - l [xg' + 2ag]y = 0. (7.2) 

Since Eq. (7.2) should be satisfied identically in the variables x, y, y', it 
splits into two equations: 

x 2 f' + axf — a(a — 1) = 0, xg + 2ag = 0. (7.3) 

Solving the first-order linear differential equations (7.3) for the unknown 
functions f(x) and g(x), we obtain: 

f( x ) = - + Ax~ a , g(x) = Bx~ 2a , A,B = const. 
x 

Thus, we arrive at the following linear equation admitting the operator 
(7.1): 

x 2a y" + (Ax a + ax 20 ' 1 ) y' + By = 0, A, B = const. (7.4) 

Remark 7.1. When a — 0, Eq. (7.4) yields the equation (6.2) with con- 
stant coefficients. When a — 1, Eq. (7.4), upon setting A + 1 as a new 
coefficient A, coincides with Euler's equation (6.16). 
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Since Eq. (7.4) is linear homogeneous, it admits, along with (7.1), the 
operator 

d 

Now we proceed as in Section 6 and find the invariant solutions for the 
linear combination X = X\ + AX2 : 

d d 
X = x a — — h Xy— , A = const. 
ox ay 

Let q^I. The characteristic equation 

dy dx 

y x a 

of the equation X(J) = for the invariants J(x,y) yields the invariant 

A 1-a 

J = y ea -i x 

Hence, the invariant solutions are obtained by setting J = C, whence letting 
C = 1 we have 

y = Qih xl ~ a , (7.5) 

Differentiating we have: 

y > = \x- a eT^ xl ~ a , y" = [\ 2 x~ 2a - \ax- a - l ]e^ xl ~ a . 

Substituting in Eq. (7.4) and dividing by the non- vanishing factor e 1 ^ 1 ' 
we obtain the following characteristic equation for Eq. (7.4): 

A 2 + A\ + B = 0. (7.6) 

Eq. (7.6) is identical with the characteristic equation (6.6) for the equation 
(6.2) with constant coefficients. 

Remark 7.2. Eq. (7.5) contains Euler's substitution (6.5) for the equation 
(6.2) with constant coefficients as a particular case a — 0. 

Remark 7.3. Using the statement that Eq. (6.1) is mapped to the equa- 
tion z" = if and only if the function 

I(x) = g(x)- 1 -f(x)- 1 -f'(x) 
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vanishes (see Lemma 6.1), one can verify that Eq. (7.5) is equivalent by 
function to the equation z" = if and only if 

A 2 - 4B = and a = or a = 2. (7.7) 

The first equation in (7.7) means that the characteristic equation (7.6) has a 
repeated root, and hence there is only one solution of the form (7.5). Then, 
using the reasoning of Lemma 6.1 we can show that the general solution to 
Eq. (7.4) with A 2 - AB = 0, a = 2, i.e. of the equation 



y"4 (- 2 + -)y' + - 4 y = o, A 2 -4B = 0, 

x z x I ar 



is given by 



V=(*i + § )e"*, (7.8) 



where Ki,K 2 are arbitrary constants and A is the repeated root of the 
characteristic equation (7.6). For a more general statement, see Section 9. 

7.2 Second example 

Let us find the linear second-order equations (6.1) admitting the projective 
group with the generator 

Taking the second prolongation of the operator (7.9), 

Xi = x 2 — + xy— + (y- xy')— - 3xy' '— , 
ox ay ay ay 



and writing the invariance condition of Eq. (6.1), 
X 1 (y" + f(x)y' + g(x)y) 

we obtain: 



= 0, 

(6.1) 



x(xf + 2f)y' + (x 2 g> + Axg + f)y = 0. (7.10) 

Eqs. (7.10) yield: 

r, X A / X B A 

f{x) = , g(x) = — - — , A,B = const. 
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Thus, we arrive at the following linear equation admitting the operator 
(7.9): 

x 4 y" + Ax 2 y' + (B - Ax)y = 0, A,B = const. (7.11) 

Eq. (7.11) is linear homogeneous, and hence admits, along with (7.9), 
the operator 

d 

X * = V W 

Now we proceed as in Section 6 and find the invariant solutions for the 
linear combination X = X\ + AX 2 : 

d d 
X = x 2 — + (x + \)y— , A = const. 
ox ay 

Rewriting the characteristic equation of the equation X(J) = for the 
invariants J(x, y) in the form 

dy x + A 

-ax 



x 1 



we obtain the invariant 



J = — e* . 



x 

Setting J = C and letting C = 1 we obtain the following form of the 
invariant solutions: 

y = xe~* . (7-12) 

Substituting (7.12) in Eq. (7.11) we reduce the differential equation (7.11) 
to the algebraic equation 

X 2 + AX + B = 

which is identical with the characteristic equation (6.6) for the equation 
(6.2) with constant coefficients. 

Example 7.1. Solve the equation 

u 2 

y" + — y = 0, u = const. (7. 13) 



This is an equation of the form (7.11) with A = 0, B = u 2 . The alge- 
braic equation (6.6) yields Ai = —iu, A 2 = iou, and hence we have two 
independent invariant solutions (7.12): 

y 1 = xe 1 *, y 2 = xe~ 1 ^. 
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Taking their real and imaginary parts, just like in the case of constant coef- 
ficient equations, we obtain the following fundamental system of solutions: 



Hi = x cos 



x 



y 2 — x sin 



(7.14) 



Hence, the general solution to Eq. (7.12) is given by 



y = x 



Ci cos (^j + C 2 sin {^j 



We can also solve the non-homogeneous equation 



UJ 



y" + —y = F(x), 



e.g. by the method of variation of parameters, and obtain 



y = x 



x 

+ - 

UJ 



(7.15) 



dcos {^j +C 2 sin {^j (7.16) 
cos ^— j J xF(x) sin ^— j dx — sin ^— j J xF(x) cos ^— j dx 



8 General result for second-order equations 
8.1 Main statements 

Note that the operators given in Eqs. (6.3), (6.17), (7.1), (7.9) are particular 
cases of the generator 



X 1 ^Hx)^ + a(x)y^ 



ID 



of the general equivalence group of all linear ordinary differential equations. 
We will find now all linear second-order equations (6.1) admitting the op- 
erator (8.1) with any fixed functions (f>(x) and a(x). 
Taking the second prolongation of the operator (8.1), 

* = *j£ +ff v%j+ Wy+(-<P')y'} [*"yH^'-r)y'+(«-W] ^ , 

and writing the invariance condition of Eq. (6.1), 



X 1 (y" + f(x)y' + g(x)y) 
23 



(6.1) 



o, 



we obtain: 

(0/' + f<t>' + 2a' - <t>")y' + {<t>g' + 2<p'g + fa' + a")y = 0. 
It follows: 

0/' + /0' + 2a'-0" = O, 
^ + 2<t>'g + /a' + a" = 0. 

The first equation (8.2) is written 

(0/)' = (0' - 2a)' 

and yields: 

f(x) = -\A + <p'-2a\, A = const. 



Substituting this in the second equation (8.2), we obtain the following 
homogeneous linear first-order equation for determining g(x) : 

( f )g ' + 2 ( t>'g = -a"-^-[A + ( f>'-2a]. 

The homogeneous equation 

4>{x)g' + 2<P'(x)g = 
with a given function <f>(x) yields 

C 

^ (j) 2 (x) 

By variation of the parameter C, we set 

u(x) 



9 = 



<t> 2 {x) ' 

substitute it in Eq. (8.4) and obtain: 

u' = -Aa' + 2aa' - 0V - <f>a" = -(Aa)' + (a 2 )' - (0a')', 

whence 

u = B — Aa + a 2 — <j>a', B = const. 

Therefore, 

Thus, we have arrived at the following result. 
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Theorem 8.1. The homogeneous linear second-order equations (6.1) ad- 
mitting the operator (8.1) with any given functions = (j>(x) and a = a(x) 
have the form 

0V+ (A + <t>' -2a)<t>y'+ (B - Aa + a 2 - <j)a')y = 0. (8.6) 
Now we use the homogeneity of Eq. (8.6) characterized by the generator 

d 

Namely, we look for the invariant solutions with respect to the linear com- 
bination X = X 1 + \X 2 : 

d d 
X = (j)(x)— + (a(x) + \)y— , A = const., 

and arrive at the following statement reducing the problem of integration of 
the differential equation (8.6) to solution of the quadratic equation, namely, 
the characteristic equation as in the case of equations with constant coeffi- 
cients. 



Theorem 8.2. Eq. (8.6) has the invariant solutions of the form 

.,. <: ^r dx , (8.7) 
where A satisfies the characteristic equation 



A 2 + A\ + B = 0. 



Proof. We solve the equation X(J) = for the invariants J(x,y), i.e. 
integrate the equation 

dy _ a(x) + A 
y <f>{x) 

and obtain the invariant 

_ r <t(x) + \ dx 

J = ye J ^ 

Setting J = C and letting C = 1 we obtain Eq. (8.7) for the invariant 
solutions. Thus, we have: 

y = e i^r dx ) y> = —L— e J^xT dx (8.9) 



y 



"=7a [(^ + A) 2 -( ( 7 + A)0 / + C t / ] e f 



4>(x) ax 



Substituting (8.9) in Eq. (8.6) we obtain Eq. (8.8), thus completing the 
proof. 
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8.2 Distinct roots of the characteristic equation 



It is manifest that if the characteristic equation (8.8) has distinct real roots 
Ai 7^ A 2 , the general solution to Eq. (8.6) is given by 



y(x) = K^ 1 ^^^ + K 2 e J ^K^ ax , K u K 2 = const. (8.10) 

In the case of complex roots, Ai = 7 + i9, A 2 = 7 — i9, the general 
solution to Eq. (8.6) is given by 



o-(x)+A 1 



dx 



J 



dx 



y(x) 



K\ cos 9 



dx A ^ . ( „ f dx 
+ K 2 sin 9 



4>{x) 



4>{x) 



dx 



(8.11) 



8.3 The case of repeated roots 



Theorem 8.3. If the characteristic equation (8.8) has equal roots Ai = \ 2 , 
the general solution to Eq. (8.6) is given by 



y 



Ki + K 2 



dx 
(f)(x) 



a{x) J r X 1 



e J *(-) , K l ,K 2 = const. 



(8.12) 



Proof. Let us new variables t and z defined by the linear first-order equa- 
tions 

dt dt dt 

X 1 (t)^ ( f>(x)— + a(x)y— = l, X 2 (t)=y— = (8.13) 

and 

X 1 (z) = <Kx)^ + *(x)y^ = 0, X 2 (z)=y^ = z, (8.14) 
respectively. Eqs. (8.13) are easily solved and yield 

t= I —■ (8.15) 



4>{x) 

Integration of the second equation (8.14) with respect to y gives 

z = v(x)y. 

Substituting this in the first equation (8.14) we obtain 



ix) — — h a(x) v = 0, whence v — e J ^> . 
dx 



tr(x) 
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Thus, 

z = ye -fj$ dx . (8.16) 

The passage to the new variables (8.15), (8.16) converts the operator 
X\ given by (8.1) to the translation generator without changing the form 
of the dilation generator X 2 . In other words, upon introducing the new 
independent and dependent variables t and z given by (8.15) and (8.16), 
respectively, we arrive at the operators (6.3). Hence, in the new variables, 
Eq. (8.16) becomes an equation with constant coefficients. Invoking that 
the equations (8.6) and (6.2) have Eq. (8.8) as their common characteristic 
equation, we use Lemma 6.1 and write 

z = (K 1 + K 2 t)e Xlt . 

Substituting this in Eq. (8.16) and replacing t and z by their expressions 
(8.15) and (8.16), respectively, and solving for y, we finally arrive at Eq. 
(8.12). 

Remark 8.1. We can easily solve the non-homogeneous equation Eq. (8.6): 

0V+ (A + (/>' -2a)<f>y'+ (B - Aa + a 2 - (f)a)y = F(x). (8.17) 

Namely, we rewrite Eq. (8.6) in the form 

y" + a(x) y' + b(x) y = P(x) 

and employ the representation of the general solution (see, e.g. [2], Section 
3.3.4) 

y = K lVl (x) + K 2V2 (x) - Vl (x) f y^^ldx + y 2 (x) f dx 

(8.18) 

furnished by the method of variation of parameters. Here 

W(x) = yi(x)i/ 2 (x) - y 2 (x)y[(x) 

is the Wronskian of a fundamental system of solutions yi(x), y 2 (x) for the 
homogeneous equation 

y" + a(x) y' + b(x) y = 0. 
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9 Examples to Section 8 

Euler's substitution (6.5) as well as the solutions (6.19), (7.5) and (7.12) 
are encapsulated in Eq. (8.7). We will consider now these and several other 
examples. 

Example 9.1. Let us take <p(x) = 1, a(x) = 0. Then Eqs. (8.6), (8.7) and 
(8.12) coincide with Eqs. (6.2), (6.5) and (6.9), respectively. Eq. (8.11) 
becomes (6.8). 

Example 9.2. Let us take <f>(x) = x, a(x) = 0. Then Eq. (8.6) becomes 
Euler's equation written in the form (6.21), Eq. (8.7) yields Eq. (6.19) for 
invariant solutions, whereas Eq. (8.12) provides the general solution 

y(x) = (K 1 + K 2 lux) x Xl K 1 , K 2 = const., (9.1) 

to Euler's equation (6.21) whose characteristic equation (6.6) has equal 
roots. Eq. (8.11) leads to the following solution for complex roots Ai = 
7 + iO, A 2 = 7 - i9 : 

y{x) = [tfi cos(# In x) + K 2 sin(# In x)] x 1 . (9.2) 

Example 9.3. Let us take <f>(x) = x a , a(x) = 0. Then Eqs. (8.6) and (8.7) 
coincide with Eqs. (7.4) and (7.5), respectively, whereas Eq. (8.12) provides 
the general solution 

y (x) = (Id + Kix^e^* 1 "*, K U K 2 = const., (9.3) 

to Eq. (7.4) whose characteristic equation (7.6) has equal roots Ai = \ 2 . 
Eq. (9.3) extends the solution (7.8) to all equations (7.4) with with the 
coefficients A, B, C satisfying the condition (6.12) of equal roots for the 
characteristic equation (7.6). 

Example 9.4. Let us take <p(x) = 1 + x 2 , a(x) = x. Then Eq. (8.6) 
becomes 

(l + x 2 ) 2 y"+ (l+x 2 )Ay' + (B-Ax-l)y = 0. (9.4) 
Working out the integral in Eq. (8.7), 

f a(x) + A , f x f A r 

/ — — - — ax = ax + ax = Inyl + x z + A arctan x, 

J <p(x) J 1 + X 2 J 1 + X 2 

we obtain the following expression for the invariant solutions: 

y = v^T^e Aarctana; , (9.5) 
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where A satisfies the characteristic equation (8.8): 

A 2 + A\ + B = 0. (9.6) 

If the characteristic equation (9.6) has distinct real roots, Ai ^ A 2 , the 
general solution to Eq. (9.4) is given by 



y(x) = VTTa^ [Ki e Al arctana; + K 2 e Aa arctanx ] . (9.7) 

In the case of complex roots, Ai = 7 + i9, \ 2 = 7 — i0, the general 
solution to Eq. (9.4) is given by 

y(x) = [Xi cos (0 arctan x) + K 2 sin (0 arctan a;)] y/l + x 2 e A arctana; . (9.8) 

Finally, if the characteristic equation (9.6) has equal roots Ai = A2, the 
general solution to Eq. (9.4) is given by 

y {x) = (Kx + K 2 arctan^Vl + x 2 e A arctan:c . (9.9) 

Example 9.5. Consider Eq. (9.4) with A = 0, B — cu 2 . Then, according 
to Example 9.4, the equation 

(l + x 2 ) 2 y"+(co 2 -l)y = (9.10) 

has the following general solution: 



y(x) = [Ki cos {uj arctanx) + K 2 sin (u arctanx)] y/l + x 2 . (9.11) 
Example 9.6. Let us solve the non-homogeneous equation 

(1 + x 2 ) 2 y" + (uj 2 -l)y = F(x). (9.12) 
Example 9.5 provides the fundamental system of solutions 



yi = Vl + x 2 cos (u arctan x) , y 2 = Vl + x 2 sin (uj arctan x) 
for the homogeneous equation (9.10). We have: 

y[ = - : = [x cos (u arctan x) — wsin (uj arctan x)\ , 



x 



y' 2 = —^=^= [x sin (uj arctan x) + uj cos (uj arctan x)\ 
Vl + x 2 

29 



Hence the Wronskian is W [2/1,2/2] = 2/12/2 ~ UiUi = u - Now we rewrite Eq. 
(9.12), in accordance with Remark 8.1, in the form 

cu 2 - 1 F(x) 

« + WT^r^jrUf (9 - 13) 

use Eq. (8.18) and obtain the following general solution to Eq. (9.12): 

y(x) = y/l + x 2 K 1 cos {ijj arctan a; ) + K 2 sin [u arctana;) 

1 f Fix) 

cos {oj arctan x) I - , sin (a; arctana:) dx (9.14) 

u> J (1 + x 2 ) 6 ' 2 

1 • / s f F ( x ) , s , 1 

H — sin w arctan x) / , cos fa; arctan x) dx . 

u> J (1 + x 2 ) 6 ' 2 J 

10 Third-order equations 

The previous results can be extended to higher-order linear ordinary differ- 
ential equations. I will discuss here the third-order equations 

2/'" + f(x)y" + g(x)y> + h(x)y = 0. (10.1) 

Theorem 10.1. The homogeneous linear third-order equations (10.1) ad- 
mitting the operator (8.1) with any given = <f>(x) and a = a(x) have the 
form 

0V'+^ + 3(0'-a)]0V (10.2) 
+ [B + A<j>' - 2Aa + <P<P" + (0') 2 - 3(0(7)' + 3a 2 ] y' 
+ [C - Bo + Ax 2 - A(f)o' -a 3 - 0V - <fx/)'a' + 30(7(7'] y = 0. 

Proof. We take the third prolongation of the operator (8.1): 

+ [ (7 "y + (2(7'-0")2/' + ((7-20')y"]A 
+ [a"'y + (3a" - 0>' + 3(a' - cf>")y" + (a - 3<P')y'"} ° 



dy'" ' 
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and write the invariance condition of Eq. (10.1): 
X 1 (y'" + f(x)y" + g(x)y' + h(x)y) 



= 0. 

(10.1) 



We annul the coefficients for y", y' and y of the left-hand side of the above 
equation and split it into the following three equations: 

4>f + 07 + 3K - 0") = 0, (10.3) 
<j>g' + 2<j)'g + (2a 1 - <f)")f - 0'" + 3a" = 0, (10.4) 
<f)ti + ty'h + a"f + a'g + a'" = 0. (10.5) 

Eq. (10.3) is written 

(0/)' = 3(0' - a)' 

and yields: 

/ = \ \A + 3((f)' -a)], A = const. (10.6) 

We substitute (10.6) in Eq. (10.4), integrate the resulting non-homogeneous 
linear first-order equation for g and obtain: 

g = -i [B + A4>' - 2Aa + #" + (0') 2 - 3(<H' + 3^ 2 ] , (10.7) 

where B is an arbitrary constant. Now we substitute (10.6), (10.7) in Eq. 
(10.5), integrate the resulting first-order equation for h and obtain: 

h = \C - Ba + Aa 2 - A(j)a' -a 3 - <j) 2 a" - 00V + 3<f)aa'] , (10.8) 

where C is an arbitrary constant. Finally, substituting (10.6), (10.7) and 
(10.8) in Eq. (10.1), we arrive at Eq. (10.2). 



Theorem 10.2. Eq. (10.2) has the invariant solutions of the form (8.7), 

y = e f^r dx , (8.7) 
where A satisfies the characteristic equation 

A 3 + AX 2 + BX + C = 0. (10.9) 

Proof. Adding to Eqs. (8.9) the expression for the third derivative y'" and 
substituting in Eq. (10.2) we obtain Eq. (10.9). 
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Chapter 3 



Using connection between 
Riccati and second-order linear 
equations 

11 Introduction 

Recall that the Riccati equation 



y' = P(x) + Q(x)y + R{x)y\ R(x) + 0, 



(11.1) 



is mapped by the substitution 



1 v! 



(11.2) 



R[x) u 



to the linear second-order equation 



u" + f(x)u + g(x)u = 



(11.3) 



with the coefficients 



Q(x) + 



R\x) 



g{x)=P(x)R(x). 



(11.4) 



Indeed, we have: 




and 



R u R u l 
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Substituting these expressions in Eq. (11.1) and multiplying by —Ru we 
obtain the equation 

R' 

u" - — v! = Qv! - PRu, 
R 

i.e. Eq. (11.3) with the coefficients (11.4). 

12 From Riccati to second-order equations 

12.1 Application to Equation (2.9) 

Applying Eqs. (11.3)-(11.4) to Eq. (2.9), 

y' = Q(x)y + R(x)y 2 , (2.9) 
we obtain the following second-order linear equation: 

U " = ( g + f ) U '- (12 ' 1} 

The integration yields: 

lnu' = J (q + ^jdx + \nd = J Qdx + \nR + \nd- 

Hence 

u' = C l R(x)^ Q{x)dx (12.2) 

and 

u = d J R(x)e f Q{x)dx dx + d- (12.3) 

Substituting (12.2) and (12.3) in Eq. (11.2) and denoting C = d/d, we 
arrive at the solution (3.6) to Eq. (2.9): 

e fQ(x)dx 

y = ~C + jR(x)eI^ dx dx' (3 ' 6) 

12.2 Application to Equation (2.10) 

The following examples clarify how to use the linearizable equations (2.10), 
y' = P(x) + Q(x)y + k[Q(x) - kP(x)]y 2 , k = const., (2.10) 
for integrating the corresponding second-order linear equations (11.3). 
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Example 12.1. If we apply Eqs. (11.3)-(11.4) to Eq. (4.8) from Example 
4.1, 

y' = x + 2xy + xy 2 , 
we obtain the following second-order linear equation: 

u" - ^2x + u' + x 2 u = 0. (12.4) 

Let us integrate this equation. Writing Eq. (11.2) in the form 

— = -R(x)y, 



u 



substituting here R(x) = x and the expression (4.9) for y, we get 

v/_ _ 1 + C + x 2 
u ~ ~ X l-C-x 2 ' 

Writing this equation in the form 

dlnu 2x 

x + 



dx x 2 + C - 1 

and integrating we obtain the following general solution to Eq. (12.4): 

u = K(x 2 + C - l)e x2/2 , C, K = const. (12.5) 

Example 12.2. If we apply Eqs. (11.3)-(1 1.4) to Eq. (4.10) from Example 
4.2, 

y ' = x 2 + (x + x 2 ) y + - (2x + x 2 ) y 2 , 
we obtain the following second-order linear equation: 

u" - (l + x)(x+ ^—^)u' + -x 2 (2x + x 2 )u = 0. (12.6) 

\ 2x + x 2 J 4 v ' 

Let us integrate this equation. Writing Eq. (11.2) in the form 

u' 

— = -R(x)y, 
u 

substituting here 

R(x) = ^(2x + x 2 ) 
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and the expression (4.11) for y, we get 



v! 





u 



l-(C + f{x + x 2 )e~* 2 / 2 dx) e- 2 / 2 



The integration yields In \u\ = In \K\ + (f>(x), where 





dx. 



1 - (C + f(x + x 2 )e~ x2 / 2 dx) e^/2 



Hence, the general solution to Eq. (12.6) has the form 



(12.7) 



where (f>(x) is given above and K is an arbitrary constant. 

Applying Eqs. (11.3)-(11.4) to Eq. (2.10) and using the solution proce- 
dure for Eq. (2.10) described in Section 4, we obtain the following general 
result. 

Theorem 12.1. The general solution of the second-order linear equation 



with an arbitrary constant k and two arbitrary functions P(x) and Q(x) 
can be obtained by quadratures. 

13 From second-order to Riccati equations 

It is manifest from Eqs. (11.4) that two coefficients f(x) and g(x) of a given 
second-order equation (11.3) do not uniquely determine three coefficients 
P(x),Q(x), R(x) of the corresponding Riccati equation (11.1). Namely, if 
we know solutions of an equation (11.3), we can solve by using the formula 
(11.2) an infinite set of the Riccati equations 



Q(x)x + 



Q'{x) - kP'(x) 
Q(x) - kP(x) 



u' + k [P(x)Q(x) - kP 2 {x)] u = (12.8) 



y' = R(x)y 2 - f(x) + 



R'(x)~ 



y + 




(13.1) 



R(x) 



R(x) 



with an arbitrary function R(x) ^ 0. 
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Example 13.1. Consider the following equation with constant coefficients: 

u" + u = 0. (13.2) 

Here / = 0, g — 1. Hence, the corresponding Riccati equation (13.1) has 
the form 

' = **w-T$i v+ m (13 ' 3) 

Substituting the general solution 

u — C\ cos x + C 2 sin x 
of Eq. (13.2) in (11.2) we obtain the following solution to Eq. (13.3): 

1 C\ sin x — C 2 cos x 

y = ■ 

R(x) C\ cos x + C-2 sin x 

If C 2 7^ we denote if = C 2 /Ci and writhe the solution in the form 

1 sin x — K cos x 
^ i?(x) cos x + if sin x 

or, upon dividing the numerator and denominator by cosx, 

1 tgx — if 

y = — . — , A = const. 13.4 

y i?(x) 1 + if tgx ' v ; 

If C 2 = the solution becomes 

ctgx 
y = ~RW 

which can be obtained from (13.4) by letting K — > 00. Thus, the general 
solution to the Riccati equation (13.5) is given by (13.4) where —00 < K < 
+00. 

In particular, taking in (13.5), (13.4) R(x) = e x we conclude that the 
equation 

y > = e x y 2 - y + e~ x (13.5) 

has the general solution 

= tgx ~ K e -* -00 < if < +00. (13.6) 
y 1 + if tgx ' - - v > 

Sometimes it is convenient to use a restricted correspondence between 
second-order and Riccati equations by writing Eqs. (11.3)-(11.4) corre- 
sponding to the Riccati equation (11.1) in the following form: 

R(x) u" - [R'{x) + Q(x)R(x)]u' + P(x)R 2 (x) u = 0. (13.7) 
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Example 13.2. Consider Euler's equation written in the form (6.21): 

x 2 u" + (A + l)xu' + Bu = 0, A,B = const. (13.8) 
Comparing the equations (13.7) and (13.8) we take R(x) = x 2 and obtain 

Q(x) = -^ P(x) = - 4 - 

x ar 

Thus, we have arrived at the following Riccati equation: 

!/* = *¥-—!/+ 4- ( 13 - 9 ) 

x 

We know that the solutions of Eq. (13.8) have the form u = x x . Substitution 
in (11.2) yields the following form of solutions to Eq. (13.9): 

y = -A. (13.10) 

Substituting (13.10) in Eq. (13.9) we obtain again the characteristic equa- 
tion (6.6): 

A 2 + AA + 5 = 0. (13.11) 

However, Eqs. (13.10), (13.11) provide only two particular solutions to Eq. 
(13.9). In order to find the general solution of the nonlinear equation (13.9), 
we have to construct the general solution u(x) of the linear equation (13.8) 
using the superposition principle and substitute u = u(x) in (11.2). 

Using the results of Section 8 on integrability of Eq. (8.6), 

(f) 2 y" + (A + <p' - 2a)(py' + (B - Aa + a 2 - <pa')y = 0, 
we can formulate the following general result. 
Theorem 13.1. The Riccati equation 

y' = R(x) y 2 - 

with three arbitrary functions R(x), (f>(x), a(x) and two arbitrary constants 
A, B is integrable by quadratures. 



A + <f)'-2a R' 
4> 



y + 



B-Aa + a 2 -^' 
R4 2 



[13.12) 
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14 Application to Ermakov's equation 

The above results on integration of linear equations 

u" + a(x)u' + b{x)u = (14.1) 

can be combined with Ermakov's method for solving nonlinear equations of 
the 

following form (see Editor's preface to Ermakov's paper in this volume): 

u" + a{x)u' + b{x)u = 4 e- 2 f a{x)dx , a = const. (14.2) 

Example 14.1. Using the solution (12.5) of Eq. (12.4) and applying Er- 
makov's method one can solve the following nonlinear equation: 

u" - ^2x + ^ju' + x 2 u = ax 2 e 2x2 u~ 3 . (14.3) 

Example 14.2. The nonlinear equation (14.2) associated with the inte- 
grable equation (8.6) has the form 

2 u" + ( A + 0' - 2a) u' + ( B - Au + a 2 - 0a') u = ^e 2 K 2 *-^ d \ (14.4) 

where and a are arbitrary functions of x, and A is an arbitrary constant. 
Eq. (14.4) is integrable by quadratures. 
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